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1 Introduction 

Integrable quantum field theories in domains with boundaries have attracted some at- 
tention of late, principally in the wake of a paper by Ghoshal and Zamolodchikov |]l|. 
Such theories can be specified in terms of ultraviolet data consisting of the origi- 
nal boundary conformal field theory together with a specification of the particular 
perturbations chosen, both in the bulk and/or at the boundary, or else in terms 
of infrared data, which might consist of a bulk S-matrix together with a set of re- 
flection factors encoding the scattering of each particle in the model off a single 
boundary. The reflection factors are constrained by various consistency conditions 
- unitarity, crossing-unitarity |j|], and the boundary Yang-Baxter [g] and boundary 
bootstrap MM equations - and it is natural to explore the space of solutions to these 
conditions, and then attempt to match them with particular perturbed boundary 
conformal field theories. While the first aspect has been studied by many authors, 
the second is comparatively underdeveloped, and provides much of the motivation for 
the work to be described in this paper. The strategy (already employed with much 
success in boundaryless situations) will be to study finite-volume spectra by two dif- 
ferent routes, one based on the ultraviolet data and one on the infrared data, and 
then to compare the results. The ultraviolet route will proceed via a generalisation 
of the truncated conformal space approach (TCSA) Q to boundary situations. This 
seems to us to be of independent interest, and should be useful in non-integrable sit- 
uations too. For brevity, we shall refer to this method as the BTCSA, for boundary 
truncated conformal space approach. On the infrared side, the main weapon will be 
the 'BTBA', a modification of the thermodynamic Bethe ansatz (TBA) equations [g] 
adapted to boundary situations, introduced in refs. |§|,|^]. We have made a study of 
the analytic structure of these modified equations, uncovering a number of surprises 
along the way. We have also found new equations, some of them encoding the ground 
state energy in regimes where the previously-proposed equations break down, and 
others the energies of excited states. These new equations are similar in form to those 
found in |^, y] for excited states in the more traditional (boundaryless) TBA. The 
agreement between BTCSA and BTBA results offers support for both, and in addi- 
tion allows us to obtain concrete predictions about the relationship between various 
parameters appearing at short and long distances. 

Thus far, most of our work has been confined to what appears to be the simplest 
nontrivial example, namely the boundary scaling Lee- Yang model, but the methods 
advocated certainly have wider applicability. In this paper the emphasis will be on the 
structure of the BTBA equations and the matching of their solutions with the results 
that we have obtained by means of the BTCSA. A more detailed description of the 



BTCSA method itself will appear in a companion paper, currently in preparation |1C]. 



2 The model 

The ultraviolet limit of the scaling Lee- Yang model (SLYM) is the A^(2/5) minimal 
model |T^, a non-unitary conformal field theory with central charge c = —22/5. The 
field content and the operator description depends on the geometry considered, but 
far from any boundary the local fields are those of the theory on a plane, and are left 
and right Virasoro descendents of the primary fields 1, of scaling dimension 0, and 
if, of scaling dimension x^p = A^+A(p = —2/5. Both of these are scalars, and the 
one non-trivial fusion rule \s ip y~ ^p = \ + ip. The conventional normalisation of (p is 

^{z) ip{w) = \z- w\''/'> + C^^ p{w) \z - wf/'> + . . . , (2.1) 

but this results in C being purely imaginary, and many other structure constants in 
the boundary theory being imaginary as well. For this reason, we have chosen the 
non-standard normalisation 

^{z) ^{w) = -\z- w\'"'> + C% ^{w) \z - «;|2/5 + . . . , (2.2) 

where we take C^p^p^p= — C,^^ to be a positive real number. 

If the physical geometry has a boundary, then a conformally invariant bound- 
ary condition (CBC) must be assigned to each component of that boundary, and 
the possible local fields on any part of the boundary depend on the particular con- 
formal boundary condition found there. Cardy has classified the possible boundary 
conditions [|l2[ and the field content on each of these can be found by solving the con- 
sistency conditions given by Cardy and Lewellen [^]. The boundary fields fall into 
irreducible representations of a single Virasoro algebra, so all that is needed to spec- 
ify the local field content of a particular conformally invariant boundary condition is 
the set of weights of the primary boundary fields. 

It turns out that there are two conformal boundary conditions for A^(2/5), which 
by an abuse of notation we shall label by 1 and <1>. The 1 boundary has only 
one primary boundary field, which is the identity 1, while the $ boundary has 
two, the identity and a field (j) of scaling dimension x^ = —1/5. As a result, there 
are no relevant boundary perturbations of the 1 boundary, and a single relevant 
perturbation of the $ boundary, by the field (j). As with equation (|2.2| ), we choose 
the normalisation of (j) to be 

</.(x)0(y) = -|x-y|2/5 + C7j^<^(y)|x-2/|i/5 + ... , (2.3) 

with C.,, = —Cfj a positive real number. Together equations (|2.2|) and (|2.3D 

determine the bulk-boundary constant ^B^ appearing in the expansion of ip{x + iy) 
on the upper half plane with boundary condition <1> at y = 0: 

v{x + iy) = '^BU{x){2yf "' + ... . (2.4) 



The only remaining independent structure constant appears in the operator product 
expansion of two boundary-changing operators, but is not needed to reproduce the 
results in this paper. 

This sketch will suffice for a description of the various boundary scaling Lee- Yang 
models that we shall be considering. Each can be regarded as a perturbation of one 
of the boundary conformal field theories just discussed. 

First, suppose that there is no boundary at all. There is only the bulk to perturb, 
and if the perturbation is to be relevant then there is only one bulk field, namely 99, 
to perturb by. The perturbed action 

/oo roo 

dy / dxip{x,y) (2.5) 

-00 J —00 

is integrable, and, for A > 0, results in a massive scattering theory with a single 
particle type of mass M, and two-particle S-matrix [0] 

S{9) = -{1){2) , (X) = ""|;if + Z\ • (2-6) 

smh(| - ^) 

The exact relationship between M and A was found in |15|. In the conventions 



implied by (O) and (2.5), it is 



M(A) = ^''''lf(mm-iro)Y'''\.n. = p.«2944 . . .)AV- , (2.7) 

^ ^ 55/16 r(2/3)r(5/6) ^ ^ ^ ^ 

Now add a single boundary along the imaginary axis x = 0. Then, as explained in, for 
example, |jl] , the S-matrix should be supplemented with a reflection factor encoding 
how the particle bounces off the boundary. There are four 'minimal' possibilities, 
each of which satisfies all of the consistency conditions entailed by the S-matrix while 
minimising the number of poles and zeroes in the strip < Ira 9 < vrPl: 

%) = -Q)(i)(i) > %) = -(i)(i)"(i)" (2.8) 

The first two were given in |||], while the second two are related to these by multi- 



plication by the bulk S-matrix. (As observed in |1£], this procedure automatically 
generates further solutions to the consistency conditions.) In fact, an easing of the 
minimality requirement will be required in order to match all of the boundary con- 
ditions that will be encountered. For the time being we just note that the reflection 
factor 

M0) = C,) (i) (!)"'(¥)"(¥) (¥) (¥)" (2.9) 



*It might be more natural to minimise the number of poles and zeroes in the narrower strip 
< Im9 < 7r/2 . In any event, we won't be imposing either version of minimality, but rather 
checking solutions directly against BTCSA data. 



is consistent with the bulk S-matrix for any value of the parameter b, and reduces to 
R(i), R[2) aiid i?(3) for 6 = 0, —2 and 1 respectively. 

At this stage there is no way of telling which, if any, of these solutions is actually 
realised as the reflection factor of a concrete perturbed boundary conformal field 
theory. Such a theory, in the geometry currently under consideration, will have an 
action of the form 

/oo /"O /'oo 

dy I dx<p{x,y) + h dy(t)B{y), (2.10) 

where ^x(2/5)+CBC i^ the action for the conformal field theory on the semi-infinite 
plane, with a definite conformal boundary condition at x = 0, and (psiy) is one 
of the boundary fields allowed by that same boundary condition. We shall denote 
the boundary $ with a term in the action h§ds(j){s) by ^{h). It is important to 
appreciate that since the bulk-boundary coupling i?^ is non-zero, the sign of h is 
important, just as the fact that the bulk three-point coupling is non-zero means that 
the sign of A in (|2.5| ) and (2.10) is important. 



As a final step, we add a second boundary to confine the system to a strip of 
finite width. The next section outlines how this situation can be analysed numerically, 
using a modification of the TCSA technique. 

3 Boundary TCSA 

The TCSA method of Q assumes periodic boundary conditions, and has to be 
adapted to our context of a strip of width R. A conformal mapping from the strip to 
the upper half plane sends the quantisation surface to a semi-circle of radius 1, with 
the left and right boundary-perturbing operators sitting on the real axis at —1 and 
1 respectively. The Hamiltonian is then given in terms of the boundary conformal 
field theory as 



H^f3iM,R) = ^ 






1-2:l / E>\ l-a:R 



+ /iL(^] %l(-1) + /iij(^) 4>k{1) 



(3.1) 



where (/? is the bulk primary field and 0l and (p-^ are the boundary perturbations (if 
any) applied to the left and right boundaries of the strip. The scaling dimensions of 
these operators are x^p, xl and xr respectively. The subscripts a and (3 of H^p^M, R) 
are included as reminders of the information residing in the left and right boundary 
conditions and perturbations, while M is related to the bulk coupling A according 



to the relation (|2.7|) . Since we are dealing here with relevant perturbations of the 



scaling Lee- Yang model, x^ is equal to —2/5, and each of (/»l and 0r is either absent, 
or is the boundary field (f) with scaling dimension x^ = —1/5. This latter option 
only arises when perturbing $ boundaries - as mentioned in the last section, the 1 
boundary does not support any relevant boundary perturbations. 

The essence of the BTCSA is to diagonalise the Hamiltonian ( |3.1D in a finite- 
dimensional subspace of the full Hilbert space, usually obtained by discarding those 
states of conformal weight larger than some cutoff. The matrix elements of the fields 
0l(— 1), 0r(1) and (/j(exp(i^)) can be found exactly, but in general the integral over 
9 in (^]^) has been performed numerically. 

If the Hilbert space of the model on a strip with conformal boundary conditions 
a and /? is 'H(^a,i3) j and the irreducible Virasoro representation of weight hisVh, then 



using ||12[ we found 

^(1,1) = ^0 , W(l,<j.) = ^-1/5 7 ^C*,*) = Vo © Vli/5 ■ (3-2) 

The full set of structure constants and correlation functions necessary to evaluate 
( |3.1D on these spaces can then be found using |jl3|, and will be given in [p^ . 

We used the BTCSA to investigate the model with a pure boundary perturba- 
tion (A=0), and also with combined bulk / boundary perturbations. In the first 
case, perturbing by a relevant boundary field while leaving the bulk massless gives a 
renormalisation group fiow which, while leaving the bulk properties of the conformal 
field theory unchanged, moves from one conformal boundary condition to another. 
The simplest context where we can study this phenomenon within the BTCSA is the 
bulk-conformal Lee- Yang model on a strip with (1,<&) boundary conditions. Per- 
turbing the <I> boundary by the (p field should provoke a flow to another conformal 
boundary condition with fewer relevant boundary operators. The only candidate is 
1, and that is indeed what we found for positive values oi h = hn. Since the bulk is 
massless, the only length scale in the problem is provided by the boundary field h, 
and the crossover occurs as a function of the dimensionless quantity r = h^'^R. The 
best signal is provided by a plot of scaling functions, which we define in terms of the 
energy levels En{h,R) by 

En{h,R) = ^Fnih'/''R) . (3.3) 

In figure la we display the gaps (F„(r) — Fo(r)) for the (1, <I>(/i)) boundary conditions 
with log(r^'^)=log(/ii?^'^) varying from —4.9 to 2.7, using the BTCSA truncated to 
98 states. We have used a logarithmic scale for to highlight the crossover region, in 
which there is a smooth flow between the TCri ^\ and Tiri j^x spectra. As should be 
clear from the plot, the levels move from the degeneracy pattern given by X-i/5 > the 



character of the V^i/^ representation, to that given by xo • 



,1/60 



X^i/dQ) 



1 + q + q^ + q-^ + 2q'^ + 2q^ + 3q^ + 3q^ + 4q^ + 5q^ + Gq^'^ + 7 q^^ + 0{q 



.12\ 



-11/60 



Xo(9) 



A2\ 



1 + q^ + q^ + q'^ + q^ + 2q^ + 2q^ +3q^ +3q^ +4q^^ +Aq^^ +0{q 

(3.4) 
With h negative the energies become complex at large (real) R, just as happens in the 
bulk if A is negated Q. However, in this case the spectrum remains identifiable. For 
large negative h, it splits into three components. There are two complex conjugate 
sets, the real parts of which are counted by the character Xo(q)j and a third set which 
is real and which is also counted by the character Xo{q)i but which appears to have a 
different asymptotic behaviour in h and to decouple in the h ^ oo limit. To illustrate 
this point, in figure lb we plot the first 53 scaling functions against r^'^ = hR^'^, 
for —12.2 < hR^'^<7A. The solid lines indicate real eigenvalues, and the dashed 
lines the real parts of complex eigenvalues. At the moment we are unsure quite what 
the interpretation of these complex eigenvalues is, or why the different components 
should be organised into representations of the Virasoro algebra. However similar 



spectra have been found before, for example see ref. |17| 




la) The excited state scaling function gaps 
as a function of log(r^/^) for the A/=0 LY 
model on a strip with boundary conditions 
(l,$(ft)), with h = {r/Rf^^. The multi- 
plicities are in parentheses. 



lb) The first 53 scaling functions Fn{r) 
plotted against r'^/^ for the M=0 LY 
model on a strip with boundary conditions 
(l,$(/i)), with h = {r/Rf^^. The multi- 
plicities are in parentheses. 



We have also used the BTCSA to investigate massless boundary flows for other 



unitary and non-unitary minimal models, with both integrable and non-integrable 
boundary perturbations. In unitary cases the spectrum remains real, but in other 
respects the picture just outlined seems to be rather general. More work will be 
needed before the full picture is clear, but one particular point is worth stressing: 
if the bulk is massless, then a boundary perturbation cannot affect this, and hence 
any boundary perturbation, integrable or non-integrable, must flow to a conformal, 
and indeed integrable, boundary condition under the renormalisation group. This 
contrasts with the generic behaviour of a bulk perturbation, where fine-tuning is 
required if the infrared limit is to be anything other than massive. It also helps to 
explain why we were able to observe massless boundary flows in the BTCSA, despite 
the errors caused by the truncation to a finite number of levels. In analogous bulk 
situations, the TCSA makes a rather bad job of modelling bulk-massless flows, even 
in situations where the oppositely-perturbed, bulk-massive, flows are captured fairly 



well - see for example section 5.2 of [18|. 



As a second application of the BTCSA, we have investigated combined bulk / 
boundary perturbations in order to unravel the connection between the reflection 



factors ( |2.8D , (|2.9| ) and particular perturbed boundary conformal field theories. To 
this end it is only necessary to consider the (1, 1) and (1, ^{h)) boundary conditions, 
but we defer the presentation of these results until later in the paper, after the 
discussion of the BTBA method. 

Finally, we have used the BTCSA to investigate the complete partition function 
of the scaling Lee- Yang model on a cylinder of length R and circumference L, with 
bulk mass M and general boundary conditions (a, (3) at the two ends of the cylinder. 
In terms of the spectrum of the strip Hamiltonian (13.11), this partition function is 



Y, expi-LEn). (3.5) 

E„i^spec{Ha(,(M,R)) 

This partition function can also be evaluated by propagating states along the length 
of the cylinder, giving the large- i? asymptotic 

Z„^(M, R, L) ~R_oo A„^(M, L) exp{-RE^{M, L)) (3.6) 

where Eq{M, L) is the ground state energy of the model on a circle of circumference 
L. The linear part of logA^p can be extracted by setting 

\og{A^p{M,L)) =\og{g^{M,L) gp{M,L)) + B^p L , (3.7) 

For the two boundary conditions we have, the two functions g-^ (M, L) and g^(h\ (M, L) 
can be expressed in terms of dimensionless combinations as 

5l(M,L) = /i(ML), g^^f,){M,L) = f2{ML,hL^/'^). (3.8) 



The numbers /i(0) and /2(0,0) are the "universal ground state degeneracies" dis- 
cussed in [^ of the corresponding UV conformal boundary conditions, i.e. 1 and 
$(0) respectively, with values as given in table |I|. The opposite, L ^ oo, limit must 
be taken with more care, as results are different for the two cases of massless and 
massive bulk. 

At M=0 (critical bulk), for /i > the function /2(0, hL^'^) governs the crossover 



between conformal boundary conditions $(0) and 1, so that we expect |1S] 



hm <7$w(0,i.) = lirn h{Q,K) = /i(0) . (3.9) 

However, for a massive bulk, we expect that the limit L ^ oo will lead to a purely 
massive boundary theory, for which the ground state degeneracy is 1, so that in all 
cases for the Lee- Yang model, 

lim gaiM,L) = 1 . (3.10) 

L — *oo 
M>0 

The quantities that we have just been discussing arise in a limit where the formula 
(|3.5D might be expected to break down, at least when truncated to a finite number 



of levels. Nevertheless, it turned out to be possible to extract both Eq{M, L)/M and 
the product ga{M,L) gp{M,L) from the numerical approximation to the partition 
function as calculated from the spectrum of the BTCSA Hamiltonian. 

As a first example, we calculated log |£^q(M, L)/M| from the BTCSA approxi- 
mation to the partition function, and compared it with the result obtained directly 
by using the TCSA for a circle. The result should be independent of the boundary 
conditions, and that is indeed what we found. The results are given in figure 2a, in 
which we compare the direct TCSA (or TEA) evaluation of log \Eq{M, L)/M\ for the 
model on a circle (solid line) with the result obtained from the BTCSA evaluation of 
■^l$(o) '^ifh a truncation to 29 states (points), and of Z^-^ with a truncation to 26 
states (open circles). As can be seen, the values for Eq{M,L)/M obtained via the 
BTCSA are indeed independent of the boundary conditions on the strip, and agree 
with the expected answers obtained by looking directly in the other channel. 

For our second example we estimated the ground state degeneracy functions 
log{ga{M,L)gf3{M,L)) for {a, (5) taking the two values (1,1) and (1,^(0)). The 
BTCSA calculations are rather more sensitive to truncation error than in the previous 
example, and we took results from various truncations up to 98 states, and then 
extrapolated in the truncation level. The numerical results at (ML) = 0.02 are 
compared with the exact results at L = in table |l|. 

We also performed a more detailed investigation of log( g^ (M, L) g-^ (M, L) ), and 
log{g'^{M,L) g^/fi\{M,L) ). These should both flow from their UV values in table 1 
to IR values of 0; but whereas there is a single flow for the (1, 1) boundary conditions, 
log((72[(M, L) (7$(;j)(M, L) ) depends on the dimensionless parameter hM~^'^. For 



No. of states 
• : 29 
o : 26 



• « « « « « » " » 




2a) log|^^(M,L)/M| plotted against log(ML) for 
the SLYM on a circle: TCSA for the system on a 
circle (line) compared with results from the BTCSA 
with boundary conditions (1,$(0)) (points) and 
(1,1) (circles). 



The b.c.s 

(«,/?) 


The exact value of 
log(5a(M,0)<7/3(M,0)) 


The BTCSA estimate for 
log(5„(M,0f)5;3(M,2^)) 


(1,1) 


ilog((5-V5)/10) = -0.64296... 


-0.643 ± 0.001 


(1,^(0)) 


ilog((5+\/5)/10) = -0.161754... 


-0.1616 ±0.0005 



Table 1: the BTCSA estimates and exact values of log( (7q,(M, 0) (7^(M, 0) 



small \hM^^'^\ this should still leave essentially a single crossover region from (1, <&) 
conformal boundary conditions in the UV to massive boundary conditions in the IR, 
but we should expect that for large hM~^'^ there are two crossover regions - first a 
massless boundary flow at Lh^'^ ~ 1 from the (l,*^') CBCs in the UV to effective 
(1, 1) CBCs, and then a crossover from these CBCs to the IR boundary conditions 
which should agree with the single crossover for the true (1, 1) case. 

In figure 2b we give plots of these functions against log(ML) for various values of 
{hM~^'^). The two straight lines in this plot are the exact UV values from table 1. 
The lowermost curve is log( gj {M, L) gj^ (M, L) ) from BTCSA at 43 states, showing 
a smooth flow from the UV value —0.64... to the IR value 0. The remaining curves 
are \og{g^{M,L) g^(^i,){M,L)) from BTCSA at 36 states. The values of hM'^/^, 
starting with the uppermost line and ending with the lowest line, are —0.65, —0.5, 
-0.4, -0.25, -0.1, 0, 0.1, 0.25, 0.5, 1.0, and 3.0 respectively. 
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2b) log( 51 (A'f , L) 51 (Af , L) ) and log( g^ {M, L) g^^h) {M, L) ) plotted 
against against log(ML) for various values of (/iM~^/^); details given in 
the text. Also shown are the exact UV values from table 1. 



There was one problem in finding numerical values for log( gj^ (M, L) g^th) (^j ^) ) > 
and this was estimating Baf3 L. It is only for large values of ML that equation ( |3.7D 
holds, and for large values of \hM~^'^\ this is outside the region of convergence of the 
BTCSA method. We were only able to find B^^,f^-. L directly from equation ( |3.7D 
for \hM~^'^\ < 0.1; for the remaining values of hM~^'^ , we estimated -B^^/^n L by 

extrapolation of the data for \hM~^/^\ < 0.0125. 

To mark this loss of accuracy, for —0.25 < {hM~^'^) < 0.25 we have plotted 
log{gi{M,L)g,^(^h){M,L)) for -6 < log(ML) < 5 with solid lines, but for the re- 
maining values of {hM~^'^) we have only given results for log(ML) < 1, and used 
dashed lines. We estimate the error in log{ g-^iM, L) (7$(/i)(M, L) ) at ML = 1 arising 
from the extrapolation of -Bj^^^^x L to be of the order of 0.01% for (hM~^/^) = 0.3, 

rising to 20% for {hM~^'^) = 3.0. For larger values of {hM~^'^) the errors in ex- 
trapolation render our data meaningless. 

For small values of \hM~^'^\ (less than about 0.3), we expect our results to be 
quite accurate. While the quantitative results are less good for the larger values, we 
do still believe that there is good qualitative agreement. We now comment on the 
results in figure 2b. There are essentially three different behaviours shown in this 
figure. 

1. For \hM~^'^\ ~ we found a single crossover to the massive case, with only 



10 



small deviations due to the boundary field. 

2. For {hM~^i^) > 1.0, we found clear signs of two crossover regions as L varied 
from to oo: the first (for small L) a massless boundary flow from the value 
with (1,$(0)) b.c.s to that with (1,1) b.c.s, and then a second region (for 
larger L) with a crossover to zero as the bulk mass scale dominates. We were 
not able to go to large enough values of {hM~^'^) to truly separate the two 
crossover regions, but we think that figure 2b is very suggestive of the behaviour 
we propose. 

3. For some value of [hM"^'^) between —0.8 and —0.6, there is a critical value at 
which log( (7-[|^(M, L)5r$(/j)(M, L) ) diverges - this was signalled by a pole in a 
rational fit to the data for -Bf^c/j) L for small \hM~^/^\. The predicted position 

of this pole is in agreement with the critical value (/icrit-^"^ ) = —0.68529 
discussed, from a different point of view, at the end of section 6 below. 

4 Boundary TBA 

For the rest of this paper we restrict our attention to cases where the bulk is massive, 
so that the bulk S-matrix is known. In situations where boundary reflection factors 
are also known (or conjectured), TBA-like equations for the ground-state energy on 
an interval have been put forward in refs. y,|3]- The derivation of these equations 
begins with the expression given in [l| for the boundary state l-Ba) corresponding to 
the boundary condition a: 



\Ba) = exp 



deK^{e)A{-e)A{e) 



|0), (4.1) 



where Ka{6) is related to the reflection factor for the a boundary condition by 
Ka{e) = Ra{f-9), and A{9) (denoted A^O) in §) is the Faddeev-Zamolodchikov 
operator creating a single particle with rapidity 6. (As in |^, we ignore the possi- 
ble presence of additional zero-momentum particles in this state. At least in some 
regimes, this will be retrospectively justified via a comparison with BTCSA data.) 

The next step is to express the partition function Z(^p{M,R,L) of the model on 
a cylinder of length R and circumference L as 

Z^p{M,R,L)^ ^{Ba\eM-RHp{M,L))\B(,)L (4.2) 

where Hp(M,L) is the Hamiltonian for the system living on a circle of circumfer- 
ence L, with periodic boundary conditions and bulk mass M, and \Ba)^ and l-B/3)^ 
are boundary states set up as in ( |4.1| ), but on the circle rather than the infinite 
line. Expanding the expressions for these states and then making a saddle-point 
approximation (bearing in mind the quantisation conditions imposed on the mo- 
menta by the periodic boundary conditions) ultimately leads to an expression for 
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— {log Za/3{M,R,l/M))/l which becomes exact in the limit I -^ oo. But in this hmit 

- the opposite to that considered in the discussion fohowing equation ( ^^ ) - the same 
quantity is given as E^ {M,R), the sought-after ground-state energy oi Hap{M,R). 
The calculational details can be found in ref. Q; the upshot is that Eq {M,R) is 
expressed in terms of the solution e{6) to the following nonlinear integral equation: 

e(^) = 2rcoshe- log Xaf3{e)-(I)*L{9) (4.3) 

where r = MR, L{9) = log(l+e-"(^)), f*g{e) = ^ JT^dO' f{e-e')g{9') , and 

Km = KU0)Kfs{-9) , m = -i^ log S{9) , (4.4) 



with S{9) the bulk S-matrix ( p.6D , and M is the particle mass ( p.7p . The solution 
to (^^) for a given value of r (and of any boundary-related parameters hidden inside 
the labels a and f3) determines a function c{r): 

6 r°° 

c{r) = — d9r cosh 9 L{9) (4.5) 

^ J —oo 

in terms of which E^ {M,R) (which from here on we abbreviate to Eq[R), when 
there is no danger of confusion) is 

Eo{R) = ^bndry + ShulkR " ^<^^^) ' (4-6) 

where i^bndry is a possible constant contribution to Eo{R) coming from the bound- 
aries, and ifbuik is the bulk energy per unit length. We will also be working with a 
dimensionless ground state scaling function F{r). With the bulk mass now nonzero, 
the discussion of BTBA results is most convenient if this is used to set the overall 
length scale, rather than the boundary magnetic field used in the definition ( |3.3| ) of 
the last section. Thus we set 

F{r) = ^Eo{R) = ^i^bndry + ^^buik - ^c{r) . (4.7) 



The equations (4.3)-(H^) are superficially very similar to the more familiar TBA 
equations found for periodic boundary conditions, but there are some important new 
features, as we now show. The identity Xaf3i9—iTT/3)Xaf3{9+iTT/3) = Aq,^(6') (which 
follows from the boundary bootstrap equations for the scaling Lee- Yang model) is 
enough to establish that the function Y{9) = exp(e(0)) solves the same Y-system as 



found in the standard Lee- Yang TBA, namely 

Y{9-iTT/3)Y{9 + iTT/3) = l + Y{9) . (4.8) 

From there, the standard periodicity property Y{9+5i7r/3) = Y(9) follows by simple 
substitution. However the consequent r-dependence of the solutions is different, 
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because the Aa^(0) term in ([4.3D gives the function £{9) a non-trivial behaviour near 
= 0, and this persists even in the r — > hniit. The Umiting shape of e{6) is a 
pair of 'kink' regions near to 6 = iblog(l/r), separated by two plateaux from what 
might fancifully be called a 'breather' region near to = 0. (This is in contrast 
to the situation in the usual TBA where there are just the two kink regions at 
9 K, iblog(l/r), separated by a single plateau of length 21og(l/r) |0].) As explained 



in 1 20 1 , the r-dependence of the ultraviolet-limiting solutions creeps in via interactions 
between the regions where the ^-dependence of £{9) remains non-trivial, and here 
these are separated by a distance log(l/r), half the value found for periodic boundary 
conditions. As a result, the regular part of the expansion of c{r) is generically a power 
series in r^'^, rather than the r^^'^ that is found when the boundary conditions are 
periodic. 

The other, 'irregular', parts of the function c{r) can be traced to the integral 



against rcosh^ in (4^). The mechanism is as described in |5|], though the arguments 
must be generalised to incorporate the effect of the central breather region. If, in 
terms of the blocks (x) defined in equation (|2.6|), the reflection factors Ra/f3 are equal 
to HxeA (^) ' then the final result is that 

c{r) = — ^^ I y sin — 1 r ^l-^r"^ -|- 'regular terms', (4-9) 




vr 



where 'regular terms' refers to the already-discussed power series in r^'^. As will 
be justified via a specific example in the next section, the scaling function F{r) is 
expected to expand in powers of r^' ^ alone about r = 0. This will reproduced by 
the BTBA result of (|4.7|) and (WM), so long as the irregular terms in (|4.9|) cancel 



against the explicit bulk and boundary energies included in (O). This requirement 
determines the constants <fbulk and E'bndry exactly: 

^bndry = ^( J^ siu ^ J M ; £^^,^ = -J-M\ (4.10) 

The value of <?buik is the same as is found when the boundary conditions are peri- 
odic [p. This is as it should be, since <Sbuik reflects a bulk property of the model. 

Finally we remark that the driving term in (|4.3| ), 2r cosh — log Xaf3{9) , is singular 
at the poles and zeroes of Xaf3{9). Therefore a set of r-independent points at which 
Y{9) either vanishes or is infinite can be read directly off (|4.3| ), at least in the strip 
—tt/3 < lm.9 < vr/SQ Points of the first sort, at which L{9) is also singular, are 



also seen in solutions of more usual TBA equations; they were called 'type 11' in |21]. 
Those of the second sort are a new feature of the boundary TBA, and will be dubbed 
'type III'. 



^beyond this strip, the y-system relation (i.i) should be used. 
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5 The (l,$(/i)) ground state in the BTBA 

We can now discuss the physical import of these results, with particular reference to 
the strip with (1, ^{h)) boundary conditions mentioned in section 3. For this setup 



the Hamiltonian (3.1) has (/)l = 0, 2—x^ = 12/5, and 1— xr = 6/5. Therefore, the 
conformal perturbation theory expansion of F{r) does have a regular expansion in 
r^'^, and so matches the behaviour derived in the last section within the BTBA. 

Next, we should decide which reflection factors describe the (l,<I'(/i)) situation. 
We found that substituting \ap{G) = ^(i){—G)Kb{0) into the basic BTBA equation 
(U) enabled us to reproduce BTCSA data for -0.68529 < h < 0.554411, the nec- 
essary values of b lying in the range —2 < b < 2. The natural conclusion, which 
will receive further support from its consistency with results to be reported in the 
remainder of this paper, is to associate boundary conditions with reflection factors 
as follows: 

1 ^ %)(0) (5.1) 

$(/i) ^ Rt{9) (5.2) 

The precise relation between h and b will be discussed shortly, but in terms of b the 
constant part of Eq{R) at large R now follows from the general result ( [4.10 ), and is: 



^bndry = ^f\/3-l + 2sin^JM. (5.3) 

Our numerical work with this particular \ap{9) was complicated by the presence, 
for — 2 < 6 < 2, of a double type II singularity in L{6) at 0=0. (At b=—2 there 
is no singularity, and at 6=2 it is of order 4.) Therefore, in the general case e~^'^' 
is singular at 0=0, and the direct numerical integration of equation ( [4.3| ) gives very 
unsatisfactory results: we estimated the accuracy for c{r) at small r to be between 
10~^ and 10~^. To alleviate this problem we defined 

ei9) = e{e) - (^logtanh |0 + log A„/3(^) , U^) = log (tanh^ \e + Kp{0)e-^^^^'^ 

(5.4) 
where q is the order of the pole at the origin (0, 2 or 4), 

and a = {!), (5 = b for this particular case. (The motivation for these redefinitions 
came from analogous manoeuvres performed in refs. |8,21] for certain excited-state 



TBA equations for periodic boundary conditions.) Equations (|4.3| ) and ( [4. 5]) can 
then be recast in the following form: 

£{9) = 2r cosh 9 -(P*L{e) , c{r) = ^^qr + ^ d9r cosh 9L{9) (5.6) 
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These revised equations are nonsingular on the real axis, and so can be solved nu- 
merically with higher accuracy. 

The attempt to go beyond the range — 2 < 6 < 2 exposes a couple of subtleties of 
the boundary TBA. We start with the situation as the point b = —2 is approached. A 
study of the singularities of L(9) in the complex 6-plane revealed a pair of Y{6) = —1 
('type r in the language of [^) singularities on the imaginary axis, at 9 = :iz9p say. 
They are confined to the segment < |Im^p| < (b+2)7r/6 by the double type II 
singularity at the origin, and type III singularities at ±i{b+2)TT/6. As b approaches 
—2 from above, the length of this segment shrinks to zero and the points zizOp are 
forced towards the origin. In order to continue round —2 to smaller real values of b, a 
deformation of integration contours is therefore required. This is just as occurs during 



the analytic continuation (in r) of ordinary TBA equations, discussed in refs. [y, 21|. 
When the deformed contours are returned to the real 6-axis, residue terms are picked 
up p, and the singularities at ±9p become 'active', in that their positions appear 
explicitly in the analytically-continued equations. These equations are: 

e{e) = 2r coshe - log A(i),(0) + log j^^ - <P*He) , (5.7) 

6 /""^ 24r 
c(r) = ^ dOr cosh OLie) + i sinh Op . (5.8) 

As in [y, we adopt the convention that 9p has a positive imaginary part after the 
real 9 axis has been crossed, and the corresponding singularity has become active. 
Its precise value appears as a free parameter in the equations, and must be fixed by 
imposing 

e{9p) = iTT. (5.9) 

These equations describe the ground state, for all real r, whenever b is in the interval 
—4 < b < —2. As before, the redefinitions ( |5.4| ) can be used to put the equations 
into a form better suited to numerical work. 

At 6=2, a different phenomenon occurs: two (inactive) type II singularities, at 



±i(6— 2)7r/6, hit the real 9 axis. As explained in section 3.4 of ref. |21], when in 
the standard TBA the real 9 axis is crossed by type II singularities, the continued 
equations can always be recast into the form that they had before the crossing, 
although in general some of the active singularities may have to be relocated from 
their analytically-continued positions. The same is true here, and the mechanism is 
as follows. After the point 6=2, the singularities at ibi(6— 2) are active and the initial 
analytic continuation of the basic BTBA equation ( [4.3| ) is therefore 

e{9) = 2r cosh^ - logX^,^,{9) - log ^|g ~ ^|^I^|^^g| - ^*L{9) . (5.10) 



Equation ( |5.7| ) arose in just the same way, though there was a different sign for 
the extra term, because of the opposite signs of the residues for type I and type II 
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singularities. This time there is no need for an equivalent of equation (|5.9| ), since the 
relevant singularities owe their existence to the factor Xn\i,[6) and so, as mentioned 
at the end of section 4, their positions are fixed irrespective of any other details of 
the solution. Now the identity 

g(g-i(b-2)7r/6) _ 
^'^^^Sie + iib-2)n/Q) - ^'-'^^^ (^•'') 

can be used to rewrite the analytically-continued equation as 

e{9) = 2r cosh 9 - log A(i)4_b(e) - </.*L(6') . (5.12) 



As promised, this has the same form as the basic BTB A equation (|4.3| ) , which applied 
before the point b=2 was passed, the only change being that the parameter b has 
been replaced by 4—6. It is straightforward to check that the expression (|4.5| ) for 
c{r) behaves in an analogous fashion under the continuation, and so the system of 
equations has a (somewhat hidden) symmetry about 6=2. Furthermore, as described 
in a related context in ||2l|, there is no reason for the 6-dependence to be in any way 
singular at this point: it is best thought of as a 'coordinate singularity', with the 
apparently discontinuous behaviour residing solely in the BTBA equations, and not in 
the functions Y and c that they encode. If the original equations had been rewritten 
using contours shifted away from the real 6 axis for all integrations, then the same 
functions would have been recovered, but from a system for which the point 6=2 had 
no special status. 

Combining the symmetry about 6=2 with the evident symmetry of the equations 
about 6=— 3 leads to the conclusion that the ground-state BTBA is periodic in the pa- 
rameter 6 with period 10. This might be a little surprising, since from ( |2.9| ) the peri- 
odicity in 6 of the driving term Xn-jf,{6) itself is 12, but is well confirmed by our numer- 
ical results. Consider now the 'regular' part of c, a function of 6 and r with an expan- 
sion in powers of r^'^. By the result just established, each coefficient in this expansion 
must be a periodic function of 6 with period 10. An excellent numerical fit for the first 
coefficient turns out to be —5.1041654269(9) sin((6+0.5)7r/5) (remarkably, all higher 
modes allowed by the periodicity appear to be absent). This coefficient can also be 
found in terms of h, and is equal to 24r(2/5)i/2 r(6/5)V2 r(4/5)~i (vrM)-^/^ h = 
1.6794178 M-6/5/i. 

Combining the two expressions yields the relation between b and h: 

h{b) = -0.685289983(9) sin f (6+0.5)7r/5) M^/^ . (5.13) 



With h{b) given by this formula we found an excellent numerical agreement between 
the BTBA and BTCSA data at all real values of 6. A particularly striking check came 
on setting 6=— 0.5 : from (|5.13| ), this should correspond to the 'pure' $ boundary 
condition, with no boundary field, and should therefore have an expansion in powers 
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of r^^'^. Thus all odd powers of r^'^ (and not just the first) should vanish in the 
regular expansion of c{r) at this point, and within our numerical accuracy that is 
indeed what we found: 

/ M n /I I 6(V3-1)(2-V2) 2^3 2 

- 8. 18x lO^^V^/^ + 0.554031116 r^2/^ 

+ 1.14x 10"^ r^^/^ - 0.0025228r2''/^ - 2.88x 10"'^ r^ + . . . (5.14) 

6 The generalisation to excited states 

We now turn to the excited states. In principle, the analytic continuation method 
of Q can be used to derive the relevant generalisations of the basic BTBA equations. 
We have yet to complete this in detail, but knowledge of the method allowed us to 
make an educated guess as to the form that the equations would take, which we then 
verified by means of a direct comparison with BTCSA data. For the one-particle 
states on the strip, the equations turned out to have the same form as those found 
in 1^] for too-particle states on a circle. With (l,<I>(/i)) boundary conditions they 
are: 

e{e) = 2r coshe - log A(i)5(0) + log f^~l'\ + log ^^^ " 4>*L{e) , (6.1) 

6 f^ 24r - 

c(r) = — dOr cosh 9 L(0) + i (sinh^o - sinh^o) , (6.2) 

e{eo) = i2n+l)Tri , e(^o) = -(2n+l)7ri , (6.3) 

where, as before, A(i)ft(0) = K ^i^{9)K},{—9) . As in the last section, we desingularised 
these equations before studying them numerically. For real r, ^o is equal to 9q, the 
complex conjugate of 6*0, and so only one of the conditions ( |6.3[ ) needs to be imposed. 
This is the generic form of the equations for — 2 < 6 < 2. 

For levels for which the one-particle excited-state equations retain the form of 
(|6.lD - (|6.3| ) as r — > oo, a large-r asymptotic can be extracted much as in P,^]. The 
result is 

E{R) - Eo{R) = Mcosh/3o (6.4) 

where /3o = Re (^o) satisfies 

2r sinh /3o - i log i?(i) {Po)Rb{f3o) = ^nk , (6.5) 

and the value of the integer k is fixed by a combination of the quantisation condition 
( |6.3D and the sign of (Im(6'o)— vr/6) : 

k = 2n + i(l-sign(Im (6lo)-7r/6)) . (6.6) 
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3a) Energy levels E{R)/M plotted against 
r for the (l,$(ft-)) boundary condi- 
tions, with 5=-1.5, /i=0.402803M6/^ 
BTBA results (points) compared with the 
BTCSA (lines). 




3b) The same as figure 3a, but showing 
the scaling functions F{r) to exhibit the 
ultraviolet behaviour. 



Equation (|6.q ) is just the 'boundary Bethe ansatz' (BBA) quantisation condition for 
the rapidity ib/3o of a single particle bouncing between the two ends of the interval. 
Further analytic results can be obtained, but for the rest of this section we will 
instead discuss some general features that emerge as h is varied. 

For h in the range —4 < b < —2, a manoeuvre similar to that seen in the last 
section for the ground state in the same region appears to be necessary, resulting in 
the appearance of a further pair of active singularities in the system, situated on the 
imaginary axis. Even for real r, there are then two independent singularity positions 
to be fixed, making the equations harder to handle numerically. We therefore leave 
this regime to one side for the time being, and move on to the range —2<b<—l. 
In this regime, the basic one-particle equations (^ 



(6.3) hold for all real r, for all 
one-particle levels. Figures 3a and 3b compare BTBA results (points) with the first 
few levels found using the BTCSA (continuous lines), at b = -1.5, h=0A02803 M'^/^ . 
Notice that the plotted points do not cover all of the lines visible on the figure: those 
missed correspond to states containing more than one particle in the infrared, and 
are presumably described by BTBA systems with more active singularities than the 
four (at ib^o and ±^o) present in (|6.1| )~ (|6.3|) . 

As b passes —1, the lowest excited state as seen in plots of BTCSA data breaks 
away from the other excited levels and dips below the one-particle threshold. Physi- 
cally the reason for this is clear: at b=—l an extra pole in Rb{0) enters the physical 
strip, signalling the appearance of a boundary bound state. After this point the 
infrared behaviour of this level changes from that of a free particle bouncing between 
the two boundaries to that of a particle trapped near to the $(/i) boundary, a state 
with asymptotic gap Mcos((6-|-l)7r/6). All of this can be seen in the BTCSA data, 
with the value of b related to h via ( ^.13| ). Since h{b) had previously been obtained 
by a matching of data in the ultraviolet, this infrared agreement provides a nontrivial 
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check on the consistency of our results. 



b=0.8, BTCSR to 67 states 




b=1.5, BTCSA extrapolated 



4) Energy levels E{R)/M plotted against 
r for 6=0.8, /i=-0.499555 M^/^. La- 
belling as in figure 3a. 




5) Energy levels E{R)/M plotted against 
r for 6=1.5, /i=-0.65175M6/5. Labelling 
as in figure 3a, with BBA results plotted 
as ( o ) . 



Figure 4 shows the situation at 6=0.8, /i=— 0.499555 M^'^, by which stage the dip 
in the first level has become quite pronounced. The points on the graph show BTBA 
results found using ( ^ ) and (|0|)-( ]6l3D , and it will be observed that the lowest 
set of excited points stops short. (The other sets persist to r=cxD and for these the 
derivation of the BBA asymptotic ( |6.4[ ) is valid.) At r ~ 4, the basic excited-state 
BTBA ( |6.3D for the first excited level breaks down, with an initial transition similar 
to that observed in [^ for the second excited state of the T2 model. This refiects the 
fact that the particle is becoming bound to the ^{h) boundary. The BTBA equations 
become more complicated in this regime, and we postpone their detailed study to 
future work. However in some regimes preliminary predictions can be obtained, in 
the spirit of [|^, by supposing that solutions to the BBA equation (|6.5| ) can be 
continued to complex values of /3o- This amounts to approximating the problem by 
the quantum mechanics of a single particle refiecting off the two walls with amplitudes 
i?(i) and Rb- For —1 < 5 < 0, the first excited level as found from the BTCSA can 
be modelled with reasonable accuracy by a solution /3o(^) to the BBA which starts 
off real at small r, reaches zero at some (6-dependent) point r^ and then becomes 
purely imaginary, tending to i(6-|-l)7r/6 from below as r tends to infinity. However 
for 5 > the only candidate solution is satisfactory at best only for large r, which 
is why such points have been omitted from figure 4. (The relevant /3o(^) tends to 
i{b+\)TT/Q from above at large r, implying an approach to the asymptotic mass gap 
from below rather than above, but our large-r BTCSA data is not yet precise enough 
to check on this prediction.) 

Continuing to increase h, the next development occurs at 6=1, /i=— 0.554411 M^'^, 
where a second excited level starts to drop below threshold, this time tending to a 
gap of Mcos((6— l)7r/6) as i? — > cxd. This heralds the arrival of a second boundary 
bound state pole in Rb{0) onto the physical strip. Figure 5 shows the situation for 
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6=1.5. For this value of b, the predicted dip in the second excited level is rather small 
(about 0.03 M ) and is only seen at rather large values of R, where the BTCSA is at 
the limits of its useful range. A more accurate reflection of the infrared behaviour 
of this second excited level is probably provided by the analytically-continued BBA 
solution plotted on the figure as open circles. The BTBA equations ( |6.3| ) now break 
down for both the first and the second excited levels, at r ~ 4 and r « 9 respectively. 
Finally, at 6=2, /i=/icrit=— 0.68529 M^'^, the asymptotic mass gap of the lowest 
excited state hits zero. It is not possible to increase h any further without making b 
complex; if a larger value is used in the BTCSA then we find energy levels becoming 
complex at large r. This can be seen analytically by combining ( ^.13| ) and (5.3) 



to write -Ehndry ^s a function of /i: a square-root singularity is revealed at h=hcrit, 
matching perfectly behaviour that can be seen directly in the BTCSA. Presumably, 
the boundary perturbation has destabilised the bulk vacuum, a situation that will 
most probably repay further study. 

7 Other combinations of boundary conditions 

Consideration of the strip with (1, <l>(/i)) boundary conditions has sufficed to establish 
the key relation ( |5.13| ), and has also enabled us to check on the basic consistency of 
our associations (|5.lD and ( ^.2| ) of the reflection factors Rti\ and Rf, with the 1 and 
$(/i(6)) boundary conditions. However, it is natural to wonder how the two other 
possible pairs of boundary conditions can be described, and this is the topic of this 
section. 

First, to the (1, 1) situation. From the identifications made earlier, one might 



expect that ground state would be described by the basic BTBA system (4.3), with 
^ap{(^) = -f^(i)(^)-f^(i)(~^)- However, a calculation of the ultraviolet central charge 
as in yields the answer c(0) = 2/5, instead of the —22/5 expected on the basis of 
the conformal result (p.2|). Moreover, numerical comparisons of BTBA and BTCSA 
results show no agreement. To motivate the resolution of this problem, we recall from 
section 3 that the 1 boundary condition can be obtained from the $(/i) boundary 
condition by sending h to +oo. That was with the bulk massless, but if we consider an 
interval of length R with (1, ^{h)) boundary conditions and take the -\-h/M^'^ — > oo 
limit while keeping R of the same order as the bulk scale 1/M, then it should be 
possible (after suitable renormalisations) to arrive at the correct BTBA. We hope to 
say more about this elsewhere, but for now we just observe, from (5.13|), that at the 



level of the BTBA the desired continuation will be found on setting 6 = —3 -|- ib, and 
then varying 6 from zero to infinity. The starting-point for this continuation lies in 
the —4 < b < —2 regime of the (1, ^{h)) ground-state BTBA, and is hence described 
by the modified system (^)-(|5^). This observation gave us the hint that the (1, 1) 
ground state might also be described by a BTBA system in which two singularities 
are already active. We therefore returned to the numerical study of the modified 
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equations Q^dH), this time with X^f^iO) = K(i)(6l)K(i)(-6'). The BTCSA data 
for the (1,1) ground state was matched perfectly over the full range of r, with a 
behaviour that is very similar to that of the first excited-state TBA for the SLYM on 
a circle [^,|^. The equations exhibit a single transition between an infrared regime, 
where two active singularities lie on the imaginary 9 axis, and an ultraviolet regime 
where these are split and join the kink systems at ^ ~ iblog(l/Mi?). As a result, 
each of these kink systems has a type II singularity sitting on the real 6 axis, a fact 
which modifies the predicted value of the ultraviolet effective central charge c(0). 
The calculation is just as described in refs. |8|,^, and the desired result c(0) = —22/5 
is indeed recovered. 

For ($(/i), <l>(/i')), calculations for the BTCSA are more involved, since two 



irreducible Virasoro representations appear in the decomposition (|3.2D of 'H($^$). 
Work on this is still in progress, but in the meantime we have considered the 
one case where a simple prediction can be made against which to test conjectures 
for the BTBA, namely h = h' = 0. With the boundary fields set to zero, the 
regular part of c{r) should exceptionally expand in powers of r^^'^, rather than 
r^'^, just as seen for the (1,<I>(0)) boundary condition in ( ^.14 ). For this case we 



found unambiguously in favour of a BTBA of the modified type ( |5.7D -(^^), with 
Ao/3(6') = i^_o.5(^)-f^-o.5(— ^)- The fit of c{r) to the two irregular terms plus a series 
in r^' ^ for this proposal yields 

^ M n A 12(y3-l)(y2-l) 2\/3 2 

c{r)\b=b'=-o.5 = 0.4 ^-^ — :^ Lr-^r 

- 3.89X 10"^^ rl + 1.79288 rf 

- 4.48X 10"^ rf - 0.414179 r¥ - 9.54x 10"^ r^ + . . . (7.1) 



As in the earlier fit ( 5.14| ), the numerical data is consistent with the coefficients of 



all odd powers of r^'^ being zero, and this supports the idea that the (<1>(0), $(0)) 
ground state is indeed described by the modified BTBA system ( |5.7| )-( ^l9| ). However 
this conclusion must remain preliminary while a comparison with BTCSA data is 
lacking, and more work will be needed before we can make any definitive statements 
about the situation for other values of h and h' . 

8 Conclusions 

We hope to have shown that the combination of boundary truncated conformal space 
and boundary thermodynamic Bethe ansatz techniques allows a detailed analysis 
to be made of integrable boundary models. Even for a theory as simple as the 
scaling Lee- Yang model, a rich structure has been revealed. Our results support the 
contention that the integrable boundary conditions of the model are 1, with reflection 
factor R/i^{9), and the one-parameter family ^{h), with reflection factors Rb{0). It 
is worth pointing out that these latter are a particular subset of the reflection factors 
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given by Ghoshal in |23] for the fowest breather in the sine-Gordon model, considered 
at the Lee- Yang point. It is natural to suppose that a boundary variant of quantum 
group reduction is at work here, and it would be interesting to explore this aspect 
further. Perhaps even more natural would be to make a link with a reduction of 
the boundary Bullough-Dodd model, since in that case the classically integrable 
boundary conditions already lie in a couple of one-parameter families |^^ . 

From BTCSA and BTBA results, and also from a simple consideration of their 
ultraviolet limits, it is clear that the 1 boundary condition is not the same as the ^{h) 
boundary condition for any finite value of h. Nevertheless, we note that Rn) = Rb=o- 
Thus the infrared data provided by an S-matrix and reflection factor alone is not 
enough to characterise a boundary condition completely. It is still possible that 
there is a special relationship between the 1 and ^{h(b=0)) boundary conditions, 
but further numerical work will be required before we can make any definite claims 
one way or the other. 

Finally we would like to reiterate that the boundary scaling Lee- Yang model was 
studied in this paper just as a first example. The methods that we have used should 
be applicable in many other situations, and we intend to report on such matters in 
due course. 
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